home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / fractal / kaos.lha / modellib / d4symmbroken_def.c < prev    next >
Encoding:
Text File  |  1990-01-16  |  2.5 KB  |  115 lines
  1.  
  2. int d4symmbroken_init()
  3. {
  4.     title_label = "D4 Symmetry Broken Oscillator";
  5.  
  6.     mapping_on = 0;
  7.     inverse_on = 1;
  8.     fderiv_on = 0;
  9.     enable_polar = 1;
  10.     enable_period = 0;
  11.     
  12.     var_dim = 4;
  13.     param_dim = 7;
  14.     func_dim = 2;
  15.     (void) malloc_init();
  16.  
  17.     var_label[0] = "x";
  18.     var_label[1] = "y";
  19.     var_label[2] = "z";
  20.     var_label[3] = "w";
  21.     var_polar_label[0] = "r";
  22.     var_polar_label[1] = "rp";
  23.     var_polar_label[2] = "theta";
  24.     var_polar_label[3] = "thetap";
  25.     param_label[0] = "damp";
  26.     param_label[1] = "force";
  27.     param_label[2] = "freq";
  28.     param_label[3] = "split";
  29.     param_label[4] = "mode1";
  30.     param_label[5] = "mode2";
  31.     param_label[6] = "cross";
  32.     func_label[0] = "Energy";
  33.     func_label[1] = "AngMom";
  34.  
  35.     param[0] = 2;
  36.     param[1] = 0;
  37.     param[2] = 0;
  38.     param[3] = 0;
  39.     param[4] = 0;
  40.     param[5] = 0;
  41.     param[6] = 0;
  42.  
  43.     var_i[0] = 0;
  44.     var_i[1] = 0;
  45.     var_i[2] = 0;
  46.     var_i[3] = 0;
  47.     var_polar_i[0] = 0;
  48.     var_polar_i[1] = 0;
  49.     var_polar_i[2] = 0;
  50.     var_polar_i[3] = 0;
  51.  
  52.     param_min[0]= -5; param_max[0]= 5;
  53.     param_min[1]= -5; param_max[1]= 5;
  54.     param_min[2]= -5; param_max[2]= 5;
  55.     param_min[3]= -5; param_max[3]= 5;
  56.     param_min[4]= -5; param_max[4]= 5;
  57.     param_min[5]= -5; param_max[5]= 5;
  58.     param_min[6]= -5; param_max[6]= 5;
  59.  
  60.     var_min[0]= -5; var_max[0]= 5;
  61.     var_min[1]= -5; var_max[1]= 5;
  62.     var_min[2]= -5; var_max[2]= 5;
  63.     var_min[3]= -5; var_max[3]= 5;
  64.  
  65.     var_polar_min[0]= -5; var_polar_max[0]= 5;
  66.     var_polar_min[1]= -5; var_polar_max[1]= 5;
  67.     var_polar_min[2]= -pi; var_polar_max[2]= pi;
  68.     var_polar_min[3]= -5; var_polar_max[3]= 5;
  69.  
  70.     f_p = d4symmbroken_f;
  71.     func_p = d4symmbroken_func;
  72. }
  73. int d4symmbroken_f(f,index,x,p,t,dim)
  74. int index,dim;
  75. double f[],x[],p[],t;
  76. {
  77.     double da,am,sp,m1,m2,cr,rs,ss,zw1,zw2,xy1,xy2,xt,y,z,w,off;
  78.     da = p[0];
  79.     am = p[1];
  80.     off = p[2];
  81.     sp = p[3];
  82.     m1 = p[4];
  83.     m2 = p[5];
  84.     cr = p[6];
  85.     xt = x[0];
  86.     y = x[1];
  87.     z = x[2];
  88.     w = x[3];
  89.     rs = xt*xt+y*y;
  90.     ss = z*z+w*w;
  91.     zw1 = cr*(z*z-w*w);
  92.     zw2 = 2.0*cr*z*w;
  93.     xy1 = cr*(xt*xt-y*y);
  94.     xy2 = 2.0*cr*xt*y;
  95.     f[0] = (-da - zw2)*xt + ( off + am - sp - m1*rs - m2*ss + zw1)*y;
  96.     f[1] = (-da + zw2)*y + (-off + am + sp + m1*rs + m2*ss + zw1)*xt;
  97.     f[2] = (-da - xy2)*z + ( off + am + sp - m1*ss - m2*rs + xy1)*w;
  98.     f[3] = (-da + xy2)*w + (-off + am - sp + m1*ss + m2*rs + xy1)*z;
  99. }
  100. int d4symmbroken_func(f,x,p,t,dim)
  101. double f[],x[],p[],t;
  102. int dim;
  103. {
  104.         double v0sq,v1sq,v2sq,v3sq;
  105.  
  106.     /* energy and angular momentum are defined in euclidean
  107.     coordinates */
  108.     v0sq = x[0] * x[0];
  109.     v2sq = x[2] * x[2];
  110.     f[0] = 0.5 * ((x[1] * x[1] + x[3] * x[3]) - p[0] * (v0sq+v2sq)
  111.         + 0.5 * (v0sq+v2sq) * (v0sq+v2sq) - p[1] * v0sq * v2sq);
  112.     f[1] =  x[1] * x[2] - x[0] * x[3];
  113.     
  114. }
  115.